Dynamical Chern-Simons Modified Gravity I: 
Spinning Black Holes in the Slow-Rotation Approximation 
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The low-energy limit of string theory contains an anomaly-canceling correction to the Einstein- 
Hilbert action, which defines an effective theory: Chern-Simons (CS) modified gravity. The CS 
correction consists of the product of a scalar field with the Pontryagin density, where the former 
can be treated as a background field (non-dynamical formulation) or as an evolving field (dynamical 
formulation). Many solutions of general relativity persist in the modified theory; a notable exception 
is the Kerr metric, which has sparked a search for rotating black hole solutions. Here, for the first 
time, we find a solution describing a rotating black hole within the dynamical framework, and in 
the small-coupling/slow-rotation limit. The solution is axisymmetric and stationary, constituting 
a deformation of the Kerr metric with dipole scalar "hair," whose effect on geodesic motion is 
to weaken the frame-dragging effect and shift the location of the inner-most stable circular orbit 
outwards (inwards) relative to Kerr for co-rotating (counter-rotating) geodesies. We further show 
that the correction to the metric scales inversely with the fourth power of the radial distance to 
the black hole, suggesting it will escape any meaningful bounds from weak-field experiments. For 
example, using binary pulsar data we can only place an initial bound on the magnitude of the 
dynamical coupling constant of £ 1,/4 < 10 4 km. More stringent bounds will require observations of 
inherently strong-field phenomena. 

PACS numbers: 04.50.Kd,97.60.Lf,04.25.-g,04.50.Gh 



I. INTRODUCTION 

Supermassive, spinning black holes (BHs) are expected 
to be present at the center of most galaxies. The gravita- 
tional field in the exterior of such bodies plays a critical 
role in the evolution of captured compact objects and in 
the emission of gravitational waves (GWs). In General 
Relativity (GR), such a field is described by the Kerr 
metric which is a stationary and axisymmetric solu- 
tion, parameterized exclusively in terms of the mass of 
the BH and its spin angular momentum. In modified the- 
ories of gravity, the Kerr metric need not be a solution 
to the field equations. For example, in Einstein-Dilaton- 
Gauss-Bonnet gravity, slowly rotating BH solutions were 
recently found that differ from Kerr [2]. A measured de- 
viation from the Kerr metric, either from electromagnetic 
observations 0,0] or purely gravitational wave ones 
can therefore provide insight into extensions of GR, or 
lack thereof. 

One particular theory that has received remarkable at- 
tention in recent years is Chern-Simons (CS) modified 
gravity 0], which also does not have the Kerr metric as 
a solution The action for CS modified gravity is de- 
fined by the sum of the Einstein-Hilbert action and a 
new parity-violating, four-dimensional correction. Inter- 
est in the model spiked when it was found that string 
theory unavoidingly requires such a correction to re- 
main mathematically consistent 0, [T3|. In the pertur- 
bative string sector, such a correction is necessary by the 
Green-Schwarz anomaly-canceling mechanism upon four- 
dimensional compactification (see eg. [ltj [HI- In gen- 
eral, such a correction arises in the presence of Ramond- 
Ramond scalars due to duality symmetries (l2T |. 



The CS correction to the action consists of the product 
of a CS scalar field, and the Pontryagin density *RR. 
The latter is defined as the contraction of the Riemann 
tensor with its dual. The dual to the Riemann tensor 
contains contractions of the Levi-Civita tensor, which is 
odd under a parity transformation, thus potentially en- 
hancing gravitational parity-breaking. Of course, this 
fact does not imply that parity-preserving solutions are 
disallowed in CS modified gravity, since as $ — > const, 
the modified theory reduces to GR. Instead, the CS cor- 
rection introduces a means to enhance parity-violation 
through a pure curvature term, as opposed to through 
the matter sector, as more commonly happens in GR. 

Two formulations of CS modified gravity exist that are 
actually independent theories on their own right, leading 
to different observables. In the dynamical formulation, 
the CS scalar is treated as a dynamical field, endowed 
with its own stress-energy tensor and evolution equation, 
while in the non- dynamical formulation, the CS scalar is 
an a priori prescribed function, and its effective evolution 
equation reduces to a differential constraint on the space 
of allowed solutions (the so-called Pontryagin constraint, 
defined as the vanishing of the Pontryagin density) . 

Most studies to-date have concentrated on the non- 
dynamical formulation of CS modified gravity. In par- 
ticular, the CS modification has been used to propose 
an explanation to the leptogenesis problem [l3|], as well 
as to explain the flat rotation curves of galaxies [Til ]. 
A short list of such studies includes some on fermion 
interactions llal . torsion [13, UM, the far-field behav 
ior MM M EH, GWs 
limit [23j . exact solutions J24 
BH perturbation theory [2 




the slow-rotation 
and Schwarzschild 
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Non-dynamical 


Dynamical 
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KMT [23], GY [26] 




Non-Axisymmetric 


AY |17,ia| 





TABLE I: Status of investigations of rotating BH solution in 
CS modified gravity. 



ical version of CS modified gravity remains completely 
unexplored. 

A common thread in the above results concerns the 
well-posedness and arbitrariness of the non-dynamical 
formulation. Indeed, in the non-dynamical framework, 
there are no well-motivated physical reasons for partic- 
ular choices of the scalar field, which is usually chosen 
ad-hoc with the goal of simplifying the equations. Fur- 
thermore, for a given choice of CS scalar the Pontryagin 
constraint significantly restricts the class of allowed so- 
lutions, even to the point where the non-dynamical the- 
ory may be over-constrained and lack a well-posed initial 
value problem [27| . In the dynamical framework the ar- 
bitrariness is reduced or even completely eliminated, as 
though one can freely prescribe the initial conditions for 
the field, the wave-like equation it must subsequently sat- 
isfy could conceivably force the field to evolve to a unique, 
late time solution independent of the initial conditions 
for a large class of spacetimes. Whether the dynami- 
cal framework offers a well-posed initial value problem is 
presently unknown. 

Only partial results exist on rotating BH solutions in 
CS modified gravity, which we classify in Tabic |l] 1 . The 
capitalized letters that appear on this table stand for 
the authors of the papers that dealt with the respective 
investigations, where a boldface font means that astro- 
physically plausible solutions were found. As one can see 
in the table, no work has been done in the dynamical 
formulation (except for some ending remarks in (26[). In 
the non-dynamical framework, solutions can be divided 
into two groups: axisymmetric and non-axisymmetric so- 
lutions. The second group is allowed in non-dynamical 
CS modified gravity because the gradient of the CS scalar 
selects a spacetime direction that need not be co-aligned 
with the BH spin, thus breaking axisymmetry. 

The first rotating BH solution was found by Alexan- 
der and Yunes (AY) [13, OH using a far-field approx- 
imation (where field point distance is assumed much 
larger than the BH mass) and linearly proportional 
to the asymptotic time coordinate t (the so-called canon- 
ical choice). This far- field solution is stationary though 
not axisymmetric, leading to a correction to the frame- 
dra ggin g effect, which was later used to constrain the the- 
ory [19|, |2fJ . The second rotating BH solution was found 



1 Henceforth, we use the term "BH solution" lightly, since all these 
investigations are in the perturbative regime, where one cannot 
formally speak of an event horizon, which is the defining feature 
of a BH. 



by Konno, et. al. (KMT) [23| . using a slow-rotation ap- 
proximation [28l . l29l | (where the spin angular momentum 
is assumed much smaller than the BH mass) and con- 
sidering both the canonical choice of i? as well as a new 
(KMT) choice. KMT showed that a slowly-rotating so- 
lution cannot exist with a canonical though with the 
KMT choice they found a stationary and axisymmetric 
solution with a CS correction that also leads to modified 
frame-dragging . 

The only study that has searched for exact, rotating 
BH solutions is that of Grumiller and Yunes (GY) [26| . 
who focused on stationary and axisymmetric line ele- 
ments for arbitrary GY argued against the existence of 
such highly symmetric BH solutions for arbitrarily large 
spins in the non-dynamical framework, because, although 
the modified field equations can be satisfied to linear or- 
der in the spin, the Pontryagin constraint generically can- 
not be satisfied to higher order. 

Both the far-field solution of AY and the slowly- 
rotating solution of KMT served greatly to understand 
the effect of the CS correction on physical observables, 
however these solutions are rather special. First, they 
have been developed purely with the non-dynamical 
framework in mind, which as mentioned suffers from po- 
tential well-posedness problems in addition to offering no 
physical or dynamical reasons for the particular choice 
of ■d. Second, these solutions do not carry over to the 
dynamical theory, as their choices for i? posses infinite 
energy, and thus do not describe a self-consistent pertur- 
bation of the GR solutions. These issues are discussed 
further in Appendix [XJ where we extend the KMT result 
and explicitly show that it is non-unique, even for the 
KMT choice of CS scalar. 

One is then led to reconsider the question of what met- 
ric represents a spinning BH in CS modified gravity with 
a well-defined CS coupling field that can be embedded 
in the dynamical framework. That is the topic of this 
paper, the first in a series attempting to shed light on 
the nature of dynamical CS modified gravity. Here we 
restrict attention to stationary, axisymmetric perturba- 
tions of the slow-rotation limit of the Kerr metric. In 
subsequent work we will relax these assumptions, which 
will most likely require numerical solution methods. 

To summarize the main results, the metric and scalar 
field describing the leading order corrections to the Kerr 
metric in dynamical CS modified gravity are: 

, 2 , 2 5 a 2 a ( 12 M 27 M 2 \ 2 n , , , 

ds~ = dsi + -- — j 1 + — h — — 7t I sin 2 Odtdfa 

K 4 /3k r 4 \ 7 r 10 r 2 / r 

5 a a cos(0) / 2M 18M 2 \ 

where ds^ is the slow-rotation limit of the Kerr metric, 
and a and (3 are coupling constants in the CS correction 
to the action (defined in Sec. HI)) , M is the BH mass and 
a ■ M is the BH angular momentum. Equation |T]) con- 
stitutes the first rotating BH solution in dynamical CS 
modified gravity, and can be thought of as a small de- 
formation of a Kerr BH with the addition of CS scalar 
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"hair" . Given that this correction is in the <f> — t sec- 
tor of the metric, there is a modification to the frame- 
dragging effect, though it is suppressed by three powers 
of the inverse radius relative to the GR solution, imply- 
ing that weak-field experiments are incapable of placing 
interesting bounds on the theory. In analogy to the post- 
Newtonian approximation (30j , the CS correction would 
correspond to a term of third post-Newtonian order rel- 
ative to the Kerr term. To this order, the location of the 
event horizon and ergosphere of the BH are unchanged. 

An outline of the rest of the paper follows. In Sec. |TT] 
we review the basics of CS modified gravity. In Sec. IIIII 
we solve the dynamical theory for a line element and CS 
coupling field describing slowly rotating BHs, valid for 
small CS coupling constants. In Sec.|lV]we describe some 
properties of the new solution, discuss related astrophys- 
ical implications, and place bounds on the strength of the 
CS coupling constants using binary pulsar data. Sec. |V] 
concludes and discusses possible future research. 

For completeness, in Appendix [5] we discuss aspects 
of approximate BHs in the non-dynamical theory: we 
describe solutions previously found, present several new 
solutions, and show that none of these solutions (new or 
old) carry over to the dynamical theory. We have also 
made an attempt to solve the CS modified equations to 
all orders in the rotation parameter with little success due 
to the complicated nature of the equations — this effort is 
discussed in Appendix [Bl 

We use the following conventions in this paper: we 
work exclusively in four spacetime dimensions with sig- 
nature (—,+,+,+) 31|, with Latin letters (a,b,...,h) 
ranging over all spacetime indices, round and square 
brackets around indices denote symmetrization and anti- 
symmetrization respectively, namely T( ah ) = ^(T a {, + Tb a ) 
and Tj a (,] = ^(T a i, — Tb a ), partial derivatives are some- 
times denoted by commas (e.g. 39/dr = d r 9 — the 
notation A( m ^ stands for a term of 0(m, n), which itself 
stands for terms of 0(e m ) or 0(£"). The Einstein sum- 
mation convention is employed unless otherwise specified, 
and we use geometrized units where G = c = 1. 



A. Basics 

Consider the action 

S = SeH + SCS + Stf + S'mat, 



where 

Scs - 
Stf - 

'S'mat 



Jv 



- / d A x*J=g~d *RR, 

Pi ii — - r ab 



(2) 

(3) 
(4) 



^x^S [g ab (V«0) (V 6 0) + 2V(0)] ,(5) 

(6) 



The first term in Eq. ([2|) is the standard Einstein- Hilbcrt 
term, while the second one is the CS correction, the third 
one is the scalar-field term and the last one describes 
additional matter sources, with £ ma t some matter La- 
grangian density. We here employ the following conven- 
tions: k" 1 = 16ttG; a and f3 are dimensional coupling 
constants, g is the determinant of the metric, V Q is the 
covariant derivative associated with the metric tensor g a b, 
and R is the Ricci scalar. The quantity *RR is the Pon- 
tryagin density defined via 



R R — RR 



*Tja cd j~}b 
n b H acd: 



where the dual Riemann-tensor is defined as 



-ft 6 = -£ -ft bef 



(T) 



(8) 



with e cde f the 4-dimensional Levi-Civita tensor. 

The CS coupling field, is a function of spacetime that 
parameterizes deformations from GR. If $ — const. CS 
modified gravity reduces identically to GR because the 
Pontryagin density is the total divergence of the CS topo- 
logical current K a 



i 



*RR, 



(9) 



where 



II. CS MODIFIED GRAVITY 



In this section, we describe the relevant aspects that 
define a complete formulation of CS modified gravity 
and establish some notation (for a detailed review see 
eg. (HI)- In the first subsection, we present a complete 
description of the CS modified action and the modified 
field equations of the theory. In the second subsection, 
we classify the modified theory into two formulations: 
dynamical versus non-dynamical CS modified gravity. 



K a = e 



abedj^n 
- bm 



dn 



y d dn 



(10) 



and T is the Christoffel connection. We can use this fact 
to rewrite 5*cs m a more standard form via integration 
by parts: 



S cs = a(tf K a )\ dv -^J^d 4 



-g (V a ti)K a . 



fill 



The first term is usually discarded since it is evaluated 
on the boundary of the manifold [33| , while the second 
term corresponds to the CS correction. 
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The modified field equations are obtained by variation 
of the action with respect to the metric and the CS cou- 
pling field: 



Rab 



-Cab 



1 

2k 



T„ 



1 



-,9abT 



(12) 
(13) 



where R a b is the Ricci tensor and □ = V Q V a is the 
D'Alembertian operator. The quantity C a b is the so- 
called C-tensor 



C ab = Vc€ cde(a VeR b) d + ^ * R d(ab)c 

where 

V a = V Q tf , Vab = V a V 6 i9. 
The total stress-energy tensor is 

T a b 



T-imat i rjvd 
1 ab + J a6> 



(14) 
(15) 
(16) 



where T^ at is the contribution from other matter sources 
(which we will set to zero here) , and T^ b is the scalar field 
contribution, given by 



(V«0) (V&0) - \gab (V a tf) (V a tf) - g ab V(ti) 



The strong equivalence principle (V a T° b 



(17) 

'mat - 0) is natu- 
rally preserved in CS modified gravity, provided the equa- 
tions of motion for d hold [Eq. (|13p ]. This is because 
when one takes the divergence of Eq. (fT2")h the first term 
on the left-hand side vanishes by the Bianchi identities, 
while the second term is proportional to the Pontryagin 
density via 



V a C Qb = ~-v b *RR. 
8 



(18) 



The equality of this contribution to the divergence of the 
scalar field stress-energy tensor leads to Eq. (flU|) . 



B. Two formulations of CS modified gravity 

The modified theory can be classified into two distinct 
formulations: dynamical and non-dynamical. The non- 
dynamical formulation is defined by setting (3 = 2 , in 
which case the field equations become 



Rab H C a b 



_ ( T^mat rr-rmat 
^ Mai ^ 



= *RR. 



(19) 



2 Traditionally, a = k when working in this formulation, but this 
is actually not necessary and we choose to leave this constant 
arbitrary. 



In vacuum, the right-hand side of the first equation iden- 
tically vanishes. The second equation, the Pontryagin 
constraint, which used to be an evolution equation for *&, 
has now become an additional differential constraint on 
the space of allowed solutions. 

In the non-dynamical framework not only does the 
Pontryagin constraint reduce the space of allowed solu- 
tions, one must also prescribe a priori the entire history 
of the CS coupling §. Once this prescription is made, the 
CS coupling field is effectively frozen and unaffected by 
any interaction. The so-called canonical choice of § is 



?9 — 



t 
(1 



-,0,0,0 



(20) 



which was chosen when the non-dynamical framework 
was first postulated Q- In hindsight, there is nothing 
truly "canonical" about this choice of CS scalar, except 
that it simplifies the modified field equations dramati- 
cally. Moreover, it is clearly not gauge- invariant, and the 
theory offers no motivation or argument why some par- 
ticular slicing of spacetime is preferred. The canonical 
choice has been found to be very restrictive, not allow- 
ing for axisymmetric rotating BH solutions [U [H, HH 
and per turbations of a certain parity of Schwarzschild 
BHs[27|. 

The dynamical formulation allows (3 to be arbitrary, 
in which case the modified field equations are given by 
Eqs. ([T2" ]) -([T7 ) . Equation (JTSJ) lifts the Pontryagin con- 
straint, and is now an evolution equation for the CS cou- 
pling field. Therefore, no restriction is imposed a priori 
on the allowed space of solutions. Instead of prescrib- 
ing the entire history of the CS coupling, one specifies 
some initial conditions for which then evolves self- 
consistently via Eq. (fTB")) . 

The dynamical and non-dynamical formulations are in- 
equivalent and independent theories, despite sharing sim- 
ilarities in the action and a few of the same solutions. 
Although one can take the limit f3 — ► at the level of 
the action to obtain the non-dynamical framework, one 
cannot expect that the same limit acting on solutions 
of the dynamical theory would return solutions of the 
non-dynamical framework. An intuitive way to see this 
is to consider ever smaller (3 parameters in the scalar 
field evolution equation. Since generically the Pontrya- 
gin density is non- vanishing, an ever smaller (3 forces an 
immense scalar field, whose bare value then couples to the 
metric via the C-tensor, leading to an equally large back- 
reaction on the geometry. Therefore, as one "freezes" the 
scalar field in the dynamical theory, one will generically 
violate its evolution equation, leading to divergences as 
exemplified by the solution we found in Eq. fl}. 

We conclude this section with a few words on the di- 
mensions of the coupling constants and scalar field. The 
choice of units of one of (a, (3, •&) constrains the units for 
the others. For example, if we wish the CS scalar field to 
have units [0] = L A , then [a] = L 2 ~ A and [(3) = L~ 2A , 
where L is a unit of length and where we have set [k] = 1. 
A more natural choice may be to require the CS scalar 
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to be dimensionless, as is done in scalar tensor theories, 
which would then require that [a] = L 2 and (3 be dimen- 
sionless 3 . Another natural choice would be to set f3 = a, 
thus putting S,g and Scs on the same footing; we would 
then have [d] = L~ 2 . Neither formulation requires that 
we choose specific units for i?, so we shall leave these ar- 
bitrary, as results in the literature have made a variety 
of different choices. 



III. ROTATING BLACK HOLES IN 
DYNAMICAL CS MODIFIED GRAVITY 

In this section we study rotating BHs in the dynamical 
formulation of the modified theory. The analytical study 
of stationary and axisymmetric line elements in this the- 
ory without the aid of any approximation scheme is a 
quixotic task. We therefore employ a couple of approxi- 
mations, and begin by defining and explaining them. We 
then proceed by to solve the modified field equations to 
second order in the perturbative expansion. 

A. Approximation Scheme 

We choose to employ two approximations schemes: 
a small- coupling approximation and a slow-rotation ap- 
proximation. In the small-coupling scheme we shall treat 
the CS modification as a small deformation of GR, which 
allows us to pose the following metric decomposition (to 
second order): 

9a b = gj$+Cg$(4) + <;*g% ) V) ) (21) 

where tr? is some background metric that satisfies the 

Einstein equations, such as the Kerr metric, while frl ($) 

and g K ab '(d) are first and second-order CS perturbations 
that depend on The book-keeping parameter £ labels 
the order of the small-coupling approximation. 

The slow-rotation scheme allows us to re-expand both 
the background and the ^-perturbations in powers of the 
Kerr rotation parameter a. The background metric and 
the metric perturbation then become 

9$} = vT+ehT + e 2 h^\ 

where the book-keeping parameter e labels the order of 
the slow-rotation expansion. We should recall here that 



Note here that we are employing geometric units where G = c = 
1, and thus, the action has units of L 2 . Had we employed natural 
units h = c = 1, then the action would be dimensionless and if 
[0] = L A then [a] = L~ A and [f3] = L~ 2A ~ 2 . 



the notation hi? stands for terms of 0(m,n), which 
in turn stands for a term of 0(e m ) and 0(£"). For ex- 
ample, in Eq. (f2"2")l , n^ " 1 is the background metric in the 

limit a = 0, while f& and h^ ^ are first and second- 
order expansions of the background metric in the spin 
parameter. 

Combining both approximation schemes we obtain a 
bivariate expansion in two independent parameter £ and 
e, which to second order is given by 

gab — 'lab +tll ab + t >"a& +£ S.' t ah +e n ab + t > a ab ) 

(23) 

When we refer to first-order terms, we shall mean terms 
of 0(1,0) or 0(0,1), while when we say second-order 
terms we refer to those of 0(2, 0), or O(0, 2) or 0(1, 1). 

What precisely are these book-keeping parameter e and 
C? The slow-rotation scheme is an expansion in the Kerr 
parameter, and thus its dimensionless expansion param- 
eter must be a/M. Therefore, a term in the equations 
multiplied by e™ is of O {{a/M ) n ). The small-coupling 
expansion parameter must depend on the ratio of CS 
coupling to the GR coupling, a/K, because such a combi- 
nation multiplies the C-tensor in Eq. (fT2|) . The definition 
of the C-tensor [Eq. {13}] clearly states that this tensor 
is proportional to gradients of the CS scalar, which it- 
self must be proportional to a//3 due to the ^-evolution 
equation [Eq. (fT3|) ]. We see then that the CS correc- 
tion to the metric will be proportional to the combina- 
tion £ = (a/ n)(a/ f3). Such a factor, however, is not 
dimensionless, and thus, it cannot be formally treated 
as a perturbation parameter. The only mass scale avail- 
able is that of the background metric, which to lead- 
ing order in the slow-rotation expansion is simply the 
BH mass. We shall then choose to normalize £ such 
that the book-keeping parameter is £ multiplies terms 
of O [a 2 /{n[3M 4 )}. 

The small-coupling scheme together with the structure 
of the modified field equations establish a well-defined 
iteration or boot-strapping scheme. From Eq. (fT2|) 
and (|13p , one can see that the source of the ^-evolution 
equation is always of lower order than the CS correction 
to the Einstein equations. Such an observation implies 
that we can independently solve the evolution equation 
for 1? first. The solution obtained for d can then be used 
in the modified field equations to find the CS correction 
to the metric. In principle this procedure can then be it- 
erated to find solutions to higher order in the expansion 
parameters. 

Let us illustrate such a boot-strapping scheme in more 
detail. To facilitate this, let us temporarily choose units 
P — k, such that is dimensionless. The small-coupling 
expansion parameter is then controlled by a only via £ = 
0[a 2 / (k 2 M 4 )]. In these units, we see then that the right- 
hand side of Eq. (fT3|) is proportional to C 1 / 2 , while the 
second term in Eq. (II2|) is proportional to £. In turn, 
this implies that $ is a Frobenius series with fractional 
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structure 



(24) 



n=0 



while the metric perturbation is a regular series in natural 
powers of C, as required in Eq. (|21|) . 

Alternatively, one could use different units for ■& that 
slightly change the order counting. For example, let 
us temporarily choose units /3 = a, which by dimen- 
sional analysis automatically implies [a] — L 4 and [9] = 
L~ 2 . The small coupling parameter then becomes £ = 
0(a/M 4 ), while the right-hand side of Eq. (JI3J is now 
proportional to £° to leading order. In these units both ■& 
and g a i, have expansions in natural powers of £, however 
the leading-order term of the former is 1/C larger than 
the latter. 

Irrespective of units, we see then that the -^-evolution 
equation is always of lower order relative to the modi- 
fied field equation, which leads to a well-defined boot- 
strapping scheme. Whether a term proportional to 
(a/M)(a//3) is of 0(1,1) or 0(1,1/2) depends on the 
choice of (3. For the purposes of order counting only, we 
shall assume that (3 <x a, though we shall leave all factors 
of (3 explicit. With this choice, aj f3 is of order unity and 
both z? and g a b have expansions in natural powers of £. 

This boot-strapping scheme is, in a sense, analogous to 
the semi-relativistic approximation [34j , where one mod- 
els extreme-mass ratio inspirals by solving the geodesic 
equations and neglecting the self-force of the particle. 
Even with this approximation, if the background is suf- 
ficiently complicated (eg. the Kerr metric), one will not 
be able to solve the field equations. The small-rotation 
scheme is then introduced such that the independent 
equations, derived from the bootstrapping scheme in the 
small-coupling approximation, can actually be solved an- 
alytically. 

Note that in the bivariate expansions we require that 
both £ and e are independently small. We shall not im- 
pose any restrictions on their relative size, except for re- 
quiring that e not be proportional to an inverse power of 
C, since this would break the above requirement. We em- 
phasize again that as is commonly done in perturbation 
theory e and C are only bookkeeping parameters and are 
not equal to a/M or £/M 4 , rather they multiply terms 
in the resultant equations which are of that order. Since 
these parameters do not carry any physical meaning, we 
can and shall set them to unity at the end of the calcu- 
lation. 



B. Slowly Rotating BH Solutions 

The slow-rotation expansion of the background met- 
ric can be formalized via the Hartle-Thorne approxima- 



tion [28|, l29( , where the line element is parameterized by 
= -f[l + h(r,d)]dt 2 
m{r, 9)} dr 2 



ds 2 



r 2 [1 



k{r, 9)} d9 2 



+ r 2 sin 2 [1 + p(r, 0)] [#- u}(r,6)dt} 2 , (25) 

where M is the BH mass in the absence of the CS cor- 
rection and / = 1 — 2M/r is the Schwarzschild factor. 
In Eq. ([25"|) . (t,r,9,(f>) are Boyer-Lindquist coordinates 
and the metric perturbations are h(r, 9), m(r, 9), k(r, 9), 
p(r, 9) and ui(r, 9). 

The metric in Eq. (|25[) has been written as in [H, [2t| , 
but clearly the metric perturbations must be expanded 
as a series both in £ and e. Keeping terms up to second 
order we have 



h(r,t 
m(r, ( 
k(r, t 
p(r,t 
uj(r,( 



= e fyl,0) + 6 C fyl.l) + e ^(2,0)) 

= £m(LO)+f( m (i,i)+e 2m (2,o), 

= e fyl.o) + 6 C fe(i,i) + e 2 fya.o)) 

= eP(LO) + fCP(i,i)+ £2 P(2,o)- 

= ew (l,0)+ f ( w (U)+ f2 W(2,0)' 



(26) 



Note that there are no terms of 0(0, 0) since these are al- 
ready included in the Schwarzschild structure of Eq. (|25j) . 
Also, we will assume that in the limit as a — > we 
uniquely recover Schwarzschild as the solution, which im- 
plies that all terms of 0(0, n) are zero. Thus, the CS cor- 
rection must be at least linear in the Kerr spin parameter 
a. From the slow-rotation limit of the Kerr metric in GR 
we can read-off the metric perturbations proportional to 
C°: 



^(i,o) — m (i,o) — fyi.o) — P(i,o) — 0, 
2Mb 

W(l, ) - —3-, 



(27) 



to first order and 



n, i2M) - —ycos-V- jj, k {20) = —a t ,rli. 

a 2 ( 2M 2 A 
P(2,o) - ^2 ( 1 + — sm 0J, « (2i0) = 0, 

2a 2 M ( 2n 2M . ,, „ 
hi 2 0-1 = . „ I cos 9 H sm 9 



(28) 



to second order. 

All fields must be expanded in the small- coupling and 
slow-rotation approximation, including the CS coupling 
field. To get a flavor of the leading-order behavior of fl 
we can return to its evolution equation [Eq. (flU|) ]. We 
see from this equation that d 2 d ~ (a/ (3) *RR, where the 
Pontryagin density vanishes identically to zeroth order 
in a/M. Thus, the leading order behavior of the CS 
scalar must be $ ~ (a/ (3) (a/M), which is clearly always 
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at least proportional to e. As discussed in Sec. IHI Al 
such a term is either of 0(1,0) ((3 = a) or 0(1,1/2) 
((3 = k), depending on the choice of /3; here we take 
the former view. Also, from the assumption that the 
Schwarzschild metric is the unique solution in the zero- 
angular momentum limit, we must have that # 0,n ) = 
for all n. The expansion for the CS scalar thus is 

= e (1 < 0) (r, 9) + e C (M) (r, 0) + e 2 (2 < o) (r, 0). (29) 

Let us now begin to apply the algorithm that we de- 
scribed earlier to solve the modified field equations, by 
focusing first on the evolution equation for the CS scalar. 
To 0(1,0) the evolution equation becomes 



(1,0) 
rr 



+ -W> l 



M 

r \ r 

72M 3 a a 



1 



(1,0) 



/3M 



COS0, 



(30) 



where we have set the potential V($) = 0. The solution 
to this partial differential equation is a linear superposi- 
tion of the homogeneous solution and a particular solu- 



tion: tfC 1 ' ) = 
is separable: 



t) 



(i.o) 

H 







(1,0) 



The homogeneous equation 



0g' o) (r,0) =#(r)#(0). 



(31) 



The partial differential equation then becomes a set of or- 
dinary differential equations for -d and i9, whose solution 
is 



tf(r) = £iF 



£ 2 F 



2' 2 
2' 2 



2M 



r 



-5/2 



-/3/2 



0(0) = S 3 P_ a/2 (cos0)+i;4Q- a /2(cose), (32) 



horizon, then must decay to a constant asymptotically 
faster than 1/r, which implies that a > 2 and (3 > 2. The 
first requirement cannot be realized for any real c\ < 1/4, 
thus forcing £i = 0, while the second requirement leads 
to ci < 0. Requiring finite total energy also implies $ 
cannot be proportional to hi(/). All these considerations 
then force us to -dv, = const. 



Now that the homogeneous solution has been found, 
we can concentrate on the particular one. One finds 



? (i,o) _ 5 a a cos(fl) 
p " 8 (3 M r 2 



2M 18M 



J2 



5r 2 



const., 



(35) 

where we can set the additional integration constant to 
zero because it does not contribute to the modified Ein- 
stein equations 4 . 

Now that the CS coupling field has been determined, 
we can search for CS corrections to the metric pertur- 
bations. Note that the stress-energy tensor for the CS 
scalar found here enters the modified field equations at 
0(2, 1), and thus it does not contribute to the metric 
perturbation. The modified Einstein equations decouple 
into two groups: one that forms a closed system of partial 
differential equations for h^ 1 ' 1 ' , to*- 1 ' 1 ), p^ 1 ' 1 ) and k^' 1 ', 
which consists of the (t,t), (r,r), (r,0), (9,8) and (</>, </>)- 
components of the modified Einstein equations, and an- 
other group that consists of a single differential equation 
for u/ 1,1 ), namely the (t, (/>)-component of the modified 
Einstein equations. 

The first group is independent of the CS coupling field, 
■d, since it arises exclusively from the Ricci tensor. One 
can verify that with Eq. (|35[) . the relevant components 
of the C-tensor for the first group vanish exactly. Since 
these metric perturbations do not constitute a CS defor- 
mation (ie. they are £ independent), we can set them to 
zero: hW = 0, = 0, jfe^ 1 - 1 ) = and = 0. 

The only remaining equation is that from the second 
group: 



where P(-) are Legendre polynomials of the first kind, 
Q(-) are Legendre polynomials of the second kind, F(-) 
are generalized hypergeometric functions, Ei are con- 
stants of integration and the coefficients 



5 = 1- VI - 4ci, /3 = 1 + VT - 4ci, 



(33) 



where C\ is the constant of integration that arises through 
separation of variables. 

Information about what reasonable constants of inte- 
gration are can be found by studying the solution i^C 1 ' ) 
in more detail. Let us first consider the far-field behavior 
of the solution, r M , in which limit 



2 sin 2 9w (1 gg l) + 3 sin 29 w^' 1 ' + 8rf sin 2 9w^' 1] 



+ 2r 2 fsm 2 9w\ 



2 (3k r 8 



x (3r 2 + 8Mr + 18M 2 ) 



(36) 



Once more, the most general solution is a linear combina- 
tion of a homogeneous solution and a particular solution. 
The particular solution is given by 



8 /3k: 



12 M 27 M" 



7 r 



10 r 2 



(37) 



d(r) ~ E 1 



M 
2r 



,-5/2 _ 



-En 



1 



M - 
2^ 



(34) 



By requiring that z9 be real, we immediately see that 
a e 5i and (3 G Si, which implies C\ < 1/4. Moreover, 
if we wish d to have finite total energy outside of the 



4 Interestingly, the behavior of a scalar field in a Kerr back- 
ground has already been studied when considering axion hair 
for Kerr [35, 36] and dyon |37|1 BHs and cosmological scenar- 
ios [3Si | , in the context of string theory. The solution found there 
is identical to the one found here in Eq. H35I I. 
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The homogeneous solution is a sum of generalized hy- 
pergeometric functions, whose argument is r/(2M) and 
depend on some separation constant c\ . For some values 
of this constant, the solution is purely real but it diverges 
at least linearly at spatial infinity, while for other values 
of this constant the solution is either complex or infinite. 
For this reason, we choose the integration constants that 
multiply these hypergeometric functions to be zero, thus 
yielding Eq. p7|) as the full solution. We note that this 
perturbation is indeed proportional to £ as expected and 
it has the correct units [u>] — since [£] = L 4 . 

The full gravitomagnetic metric perturbation to linear 
order in £ and e is 5 



2Mb 



5a£ 



12 M 27 M 2 



-3l 1 + TT + io72rl- (38) 



The above formulae constitute the first slow-rotating BH 
solution in dynamical CS modified gravity. Note that the 
perturbation is highly suppressed in the far field limit, 
decaying as r~ 6 , which suggests that its signature can 
only be observed in the strong field regime. 

As expected, the correction to the metric is a small 
^-deformation of the Kerr metric, in agreement with the 
small-coupling approximation. In fact, one could verify 
that such an approximation is self-consistent by calcu- 
lating the next order correction to This correction 
consists of both $( 2 >°) and i?' 1 ' 1 ), which can be computed 
by solving the evolution equation to next order. Doing 
so, we find 



0<U) 



+ 



25 a £a cos9 
~448/?M5 r 2 
80M 4 144M 5 



7r 4 



7r 5 



2M 
r 

112M 6 

5r 6 



18M 2 32M 3 



+ 



5r 2 
448M 7 
25r 7 



5r 3 
(39) 



which is clearly £-times smaller than $(i>°), thus render- 
ing the small-coupling approximation self-consistent. If 
we were to use this improved $ solution in the modified 
field equation, we would find a correction to the metric 
proportional to £ 2 e, which we are here neglecting. 



IV. PROPERTIES OF THE NEW SOLUTION 

We now wish to study some of the geometric properties 
of the slowly-rotating solution just found, and comment 
on astrophysical implications of this. For completeness, 



5 We leave k explicitly here to keep track of relative dimensions 
between a, (3 and n and also to keep track of hidden factors of 
16tt. 



we present the non-vanishing metric components below: 
2a 2 M 



9tt = -f 



• cos 2 9, 



2Ma 



9t<j> — 



■ sin 



5 £ a M 5 ( 12M 27M 2 \ 2 
+ 8W 4 M—{ 1 + — + -[& 2 -) Sm 6 > 

1 fl2 ( 2n l \ 



r 2 + a 2 cos 2 i 



= r 2 sin 2 9 + a 2 sin 2 9 ( 1 + sin 2 ■ 



2M 



(40) 



and recall that this metric is accurate to 0(2, 0), 0(1, 1) 
and 0(0,2). 

Initially one may wonder whether the CS correction 
can be gauged away via a coordinate transformation, 
though after some consideration it is clear this is not 
the case. In particular, one can calculate curvature in- 
variants and show that these are indeed CS corrected to 
the order considered here. Perhaps the most obvious of 
these invariants is the Pontryagin density, *RR, which 
must be proportional to and so the deviation from 
Kerr can readily be computed from (|39p . Also, as we 
shall see shortly, the location of the inner-most stable 
circular orbit is CS corrected, which is a further indica- 
tion that the CS modification is a non-trivial geometric 
perturbation of Kerr. 

The location of the ergosphere and event horizon are 
however unchanged by the CS correction to the order 
considered here. The ergosphere is the location where 
the Killing vector (d/dt) a becomes null, or equivalently 
when g tt = 0, and this component of the metric is 
unaltered by the CS-correction. The event horizon of 
the metric can be found by tracing an outgoing quasi- 
spherical lightcone backwards in time from X + . Such 
a lightcone is defined by an axisymmetric null hyper- 
surface u(t, r,9) = t — R(r, 9) — const, that satisfies 
u, a u,/3g al3 — [39] . It is straight-forward to check that, 
to this order in the expansion parameters, the lightcone 
equation does not depend on the CS correction, and thus 
will have the same solutions as Kerr; hence the horizon 
location is unchanged. 

An immediate consequence of the rapid 1 /r 4 decay of 



the perturbation to the Kerr metric (|40|) is that to this or- 
der the asymptotic structure of the solution is unaltered, 
and thus, the CS corrected solution has the same ADM 
mass and angular momentum [40(. A subtle point, how- 
ever, is that the "physical" mass and angular momentum 
need not be given by the ADM quantities in alternative 
theories of gravity. In order to determine the former, one 
needs to perform a study of the conserved Noether cur- 
rent and charges in the modified theory a la 24], which 
we shall not carry out here. Moreover, note that even 
though we colloquially called the scalar field part of the 
solution "hair" , it does not constitute a violation of the 
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no-hair theorems in that we cannot freely choose to 
this order in the perturbation ?? is uniquely given once 
we choose M and a. 

Next, we compute the total energy carried by the scalar 
field, the expression for which is 



Efl — 



1 al3 l 1 



(41) 



where £ is a t — const, hypersurface, t a = (d/dt) a , and 7 
is the determinant of the metric intrinsic to S. Integrat- 
ing the energy for the solution [Eqs. (|35|) - (|40f ] outside 
the horizon gives 



1685ttk a 2 £ 
36864 Af3' 



(42) 



which is of 0(2, 1) and thus beyond what is considered 
here. Thus, to this order in the expansion, the scalar field 
contributes zero total energy to the spacetime, which is 
consistent with both the horizon area and ADM mass 
not changing (and similar for the angular momentum), 
assuming the usual relationships between the geomet- 
ric properties and mass/angular-momentum of the BH 
holds. Note that if we require positive energy for the 
coupling constant £ must likewise be positive. 

We next turn to the geodesic structure of the space- 
time. One interesting modification is the location of 
the innermost stable circular orbit, or ISCO, for equa- 
torial orbits. This can be calculated as the circular orbit 
with minimum energy E. The pure GR result is given 

by mm 

Risco 



M 



Zi 

Z 2 



3 + Z 2 t V(3- Z 1 ){3 + Z 1 +2Z 2 ) 
l + (l-a 2 ) 1/3 [(l + a) 1/3 + (l-«) 1/3 
(3a 2 + Z 2 ) 1/2 . 



(43) 



To leading order in the expansion parameters, we find 
that in CS gravity the location of the ISCO is at 



R 



ISCO 



„„, 4V6a 7a 2 77-y/6a£ 

6Mt— ± — - — 7] 44 

3 18M 5184A/ 4 ' v ' 



where the upper signs are for co-rotating geodesies, and 
the lower signs for counter-rotating ones. As discussed 
above, if we require that the CS scalar 1? has a local stress- 
energy tensor that satisfies the usual energy conditions, 
then £ must be positive, and therefore the CS correction 
always opposes the change in the location of the ISCO 
relative to the leading order change that rotation intro- 
duces. I.e., for co-rotating geodesies, the CS correction 
has the effect of enlarging the ISCO, and vice-versa for 
counter-rotating geodesies. That the ISCO location is 
changed could offer potentially observable consequences 
in accreting BH systems, as the inner properties of the 
accretion disk are strongly affected by the location of the 
ISCO. 



The new solution also modifies the dragging of iner- 
tial frames by the rotation of the BH. This can most 
easily be gauged by looking at the angular velocity u>z 
of zero-angular-momentum observers (ZAMOs); loz — 
—gt<p/g<j><j>, or explicitly to leading order in the expansion 
parameters: 



2Ma 5 a £ M 5 



12 M 27 M 2 



(45) 

As with the effect on the ISCO location, if £ is posi- 
tive, the CS correction always reduces the magnitude of 
frame-dragging relative to Kerr for geodesies that ap- 
proach close to the horizon. 

For a more quantitative exploration of the CS cor- 
rection to the frame-dragging mechanism we can apply 
the gravitomagnetic formalism. Let us then first trans- 
form the metric to Cartesian coordinates via the standard 
transformation 



V = 



r ( 1 + ^"-j ) cos s ^ n ^' 



/ ( 1 + ' si " ' ,sm ' 

r cos 6, 



(46) 



expanded to leading order in Q(a/M). We can now define 
the trace-reversed metric perturbation 



h a b — h a b — ^Vabh, 



(47) 



where h = h a a is the trace of the metric perturbation, as 
well as the gravitomagnetic potential via A a — —hao/'i, 
where stands for the time-component. With this at 
hand, we can then define the gravitomagnetic field via 



which for the metric found here reduces to 
3 aM x z 15 a x z £ 



(48) 



ly. 



B,, 



2 R 2 R 3 

3 aM yz 

2"R2"^3 

laM 
2~R3 



16 R 8 

15 a y z £ 

16 R 8 ' 

3£f 

R 2 



B * = -0-55- U--55- +ol£> ( 49 ) 



5 a£ 
~RS 



where R 2 = x 2 + y 2 + z 2 . One can check that the GR 
terms are indeed the standard ones (see eg. [13|), being 
proportional to J = Ma, while the CS correction is pro- 
portional to J/M, since it is mass- independent. From 
this, we can see that the modification is comparable to 
a third post-Newtonian correction to the frame-dragging 
effect, ie. the CS correction is suppressed by a factor of 
1/R? relative to the GR correction, which for objects 
on circular orbits corresponds to a correction of 0(v e ), 
where v is the circular orbital velocity. 

Again, the preceding calculation indicates that dynam- 
ical CS modified gravity is not subject to meaningful 
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bounds via weak-field tests. In order to explicitly ver- 
ify this expectation, we can study some observable de- 
rived from the CS correction. Here, following [2(|, we 
calculate the variation of orbital elements of a binary, 
averaged over one period, by solving the Gaussian pertur- 
bation equations in the small eccentricity limit. Strictly 
speaking, these results can only be applied in the limit 
of extreme mass ratios, with the massive object a slowly 
rotating BH. We suspect that the modified Kerr solution 
will also provide a good description to the dynamical- 
CS induced correction to the exterior field of a rotating 
compact matter object, such as a neutron star, though of 
course a solution for that scenario would need to be found 
to verify the claim. We find that the relative average vari- 
ation of any orbital element, A, can be decomposed into 
A = Acr + Acs, where the first term is the GR expec- 
tation and the second term is the CS correction. We find 
then that only the average rate of change of the longitude 
of the line of nodes, f2, and the argument of the perigee, 
lu, are CS corrected, namely 



cs 



^GR 



25 £ 
64 MrL 



1 — - cos 2 I 



(50) 



where r sm is the semi-major axis, i is the inclination 
angle and the angled-brackets stand for the average over 
one orbital cycle. 

It then follows that a measurement of I and f2 or lu to 
an accuracy 8 can be used to test CS modified gravity, 
as done for example in [lj| [2(| ■ Let us assume that the 
measurement is in full agreement with the GR expecta- 
tion up to experimental uncertainties, and that i = tt/2. 
Then, one can constrain the coupling strength of the CS 
correction to be 



* < -Mrl 



(51) 



We can go one step further and use the data derived from 
observations of PSR J0737-3039 A/B [H[ to place the 
first bound on dynamical CS modified gravity, namely 
r sm = 4.2410 5 km, i rj 86 degrees and M = 1.476M . 
We then find 



£ 1/4 < 10 4 km, 



(52) 



where we have used an extremely conservative measure 
of the cumulative error 8 ~ 0.1 degrees per year. 

Although this is the first bound on dynamical CS mod- 
ified gravity, we wish to stress that there are several 
caveats that must be taken into account. First, as men- 
tioned above, we have used a BH solution to represent 
the exterior gravitational field of a neutron star. Such 
an identification can be shown to be accurate for weakly- 
gravitating systems in GR, however in CS modified grav- 
ity this need not be the case. Second, in order to use 
binary pulsar observations to test an alternative theory 
of gravity one would have to calculate at least two other 
observables, one of which is of course the rate of change of 



the semi-major axis due to GW emission. The generation 
of gravitational radiation in CS modified gravity has not 
yet been worked out, and thus, this aspect cannot yet be 
incorporated into (|52[) . In all then, the constraint should 
be considered as a conservative order-of-magnitude esti- 
mate on the the dynamical formulation of CS modified 
gravity. 



V. CONCLUSIONS 

In this work we have presented a first study of the 
nature of slowly rotating BHs in dynamical CS modified 
gravity, to leading order in the coupling constant. In con- 
trast to solutions in the non-dynamical framework, dis- 
cussed in Appendix [XJ no ad-hoc prescriptions need be 
made regarding the nature of the CS field Instead, we 
only demand the physically reasonable requirements that 
# respect the symmetries of the spacetime, and posses fi- 
nite, positive energy exterior to the BH event horizon. 
The resultant solution represents an inherently strong- 
field perturbation of Kerr, in that the deformation of 
the background geometry decays as (M/r) 4 . As such, 
the deformation can easily be consistent with all existing 
weak-field bounds, yet still offer the exciting possibility 
of allowing very different phenomena in strong-field sce- 
narios involving spinning BHs: compact object mergers, 
the inner edges of accretion disks, gravitational collapse, 
etc. 

Regarding the nature of the solution, we found that to 
the order of the perturbation studied here, the horizon 
and ergosphere locations are unaltered relative to Kerr, 
and similarly the ADM mass and angular momentum of 
the spacetime is unchanged. The dynamical field § can be 
thought-of as imbuing the BH with dipole scalar "hair" , 
though the structure of which is uniquely determined by 
the spin and angular momentum of the BH. The main 
observable effect of the CS correction is to alter the near- 
horizon geodesic structure of the spacetime, effectively 
weakening the frame-dragging phenomenon, and moving 
the ISCO farther out (closer in) for co-rotating (counter- 
rotating) orbits relative to Kerr. 

We further calculated the effect of the CS correction on 
the evolution of the parameters describing a binary orbit 
relative to GR. We then used these results with binary 
pulsar observations to place a bound on the strength of 
the coupling constant in the dynamical theory of < 
10 4 km. The caveats with this application of the new 
solution are that we have not yet studied whether it also 
gives an adequate description of the CS correction to the 
geometry exterior to a neutron star, that the orbit can 
be depicted as a geodesic about a slow-rotating BH, and 
that to this order the CS modification to GW emission 
in the system is negligible. 

Future work includes continued study of the nature 
of the dynamical theory, in particular BH solutions and 
compact object mergers, without imposing any slow- 
rotation or small- coupling limit approximations. This 
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will likely require numerical solution of the CS modified 
Einstein equations. One important question that such 
studies could shed light on is how does the weak-field 
nature of a well-defined metric theory of gravity restrict 
the nature of the strong field regime for astrophysically 
relevant solutions. 

Such a study is of particular importance to the budding 
field of GW astronomy. On the one hand GW detectors, 
such as LIGO (Laser Interferometer Gravitational Obser- 
vatory) [il, [3, HI] or the planned space-based detector 
LISA (Laser Interferometer Space Antenna) [13, 0, 
could be used to place restrictions on potential CS devia- 
tions from GR. On the other hand, such studies could also 
help us understand how seriously GW detection might 
be negatively impacted by search strategies built around 
the assumption that GR is the correct theory of clas- 
sical space and time; i.e. , could classes of strong-field 
GW sources be misidentified, or even missed altogether 
if there are significant deviations from GR in the strong- 
field regime? 

Investigations of this type could also help clarify how 
robustly conclusions can be made about the nature of 
strong-field sources. For example, let us assume that we 
know BH mergers "look" similar in terms of GW emission 
for a large class of metric theories that only differ from 
GR in the strong-field regime. Then, following detection 
of such a source, one can have significant confidence in a 
claim that black holes have been directly observed, even 
if one cannot claim they are Kerr black holes. 

Note added in revision: After this paper was sub- 
mitted, another work [f5(| appeared on the arXiv that in- 
dependently derived the solution presented in Sec. IIII Bl 
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APPENDIX A: ROTATING BH SOLUTIONS IN 
NON-DYNAMICAL CS MODIFIED GRAVITY 

In this appendix, we explore and generalize previ- 
ous results regarding rotating BH solutions in the non- 
dynamical framework. First, we review the slow-rotating 
solution of KMT and extend it to much more general CS 
scalar fields. Second, we find a new solution for the met- 
ric tensor with the KMT choice of CS scalar, thus proving 
that the KMT solution is not unique. We also show that 
the solutions found in the non-dynamical framework do 
not remain solutions in the dynamical theory. 



In the non-dynamical framework, KMT used the slow- 
rotation approximation to find a stationary and axisym- 
metric solution, provided the CS coupling field was cho- 
sen to be non-canonical. In particular, the KMT choice 
of 1? is 



KMT 



r cos 9 



(Al) 



where Ao is a constant. This choice of CS scalar leads to 
a solution where all metric perturbations vanish, except 
for the gt,j) — —r 2 uj components, where 



-KMT 



(r) = /+ ^ [r 2 f - 4Af 2 + 4Mr/ln(r/)] 



(A2) 

and Bi are integration constants. 

This solution was later used to propose an explana- 
tion to the flat rotation curves of galaxies [14| as follows. 
Consider the transverse, circular velocity of point parti- 
cles in this background vj,, in the geodesic approxima- 
tion and neglecting t?-backreaction. In the far- field limit, 
M/r -C 1, and using the metric found above [14| . 



h const., (A3) 

r 



" r 




where the first term is the Schwarzschild term and the 
second term is a CS correction, proposed as an explana- 
tion to the galactic rotation curves. 

The KMT solution is an interesting result that sheds 
some light on the effects of CS modified gravity on certain 
observables, though before one can claim victory over the 
rotation curves, their analysis must be considered more 
carefully. In doing so, we have found that the KMT 
solution suffers from a few drawbacks that render it un- 
physical as a true BH solution. The main problem stems 
from the fact that this solution was found initially in the 
non-dynamical formulation, which as we argued in the in- 
troduction is quite contrived, arbitrary and probably not 
well-posed. The hope would then be to embed this solu- 
tion in the dynamical framework, however such a task is 
impossible as we explain below. 

In the dynamical formulation the nature of the CS field 
i? changes in the following fundamental way: ?9 becomes 
a dynamical field governed by a scalar field Lagrangian, 
and consequently provides a new stress-energy contribu- 
tion to the CS modified Einstein equations. The KMT 
solution indeed satisfies the evolution equation for $ to 
leading order in the spin 14|. The stress-energy tensor 



associated with $kmt, however, has infinite total energy, 
and therefore the KMT metric is not a self-consistent so- 
lution to the dynamical field equations. 

In view of these problems, one could wish to general- 
ize the slow-rotation results in (23|, by considering more 
general CS scalars. A straightforward calculation reveals 
that the most general 1? for which the solution found 
in [23[ persists is 

i9 s ° n = Ao + A x r cos <j> sin + A y r sin <fi sin 9 + A z r cos 9, 

(A4) 
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where At are constants. In fact, we can rewrite this CS 
coupling field as = S a bA a x b , where x a = [l,x,y,z] 
and 5 a b is the Euclidean metric. Obviously, although 
more general, this scalar contains the same problems as 
$kmt, namely its associated stress-energy leads to infi- 
nite energy, making it incompatible with the dynamical 
framework. 

Additionally, one can also show that there is another 
solution to the modified field equations that was missed 
in [23I ] . In particular, a direct calculation reveals that 

= f(r,<f>) +rg{cf>) +rh(C 1( f>-t) + rk(6,(t>) 



to 



dr 
r 



-d r f(r,4>) 



-/>, <f>) + ~3(r) 
r r 



(A5) 



also solves the modified field equations, where /, g, h, 
j and k are arbitrary functions and C\ is another in- 
tegration constant. The CS coupling field found here 
does bypass the problem with the stress-energy tensor, 
since the arbitrary functions can be freely chosen to de- 
cay fast enough. For example, iff = g = h = k = 
and j = —3jo/r 2 , then — jo/r 2 , for some constant 
jo. Thus, although the canonical is not allowed by 
this family of solutions, Eq. (| A5|) is compatible with the 
dynamical framework. 

This better-behaved and much more general CS cou- 
pling, however, leads to a solution that differs sufficiently 
from the KMT one to lose its appeal. In other words, the 
new solution in Eq. (|A5[) cannot explain the galactic rota- 
tion curves because the logarithmic term has disappeared 
from the metric. One can show by direct calculation that 



V<f> 



M 



(A6) 



in the far-field limit M/r <C 1. This result is to be con- 
trasted with Eq. (|A3p . which leads to a constant term 
instead of a term that decays quadratically with radius. 
Due to this decaying behavior, the new solution cannot 
explain the flat rotation curves. 

One could argue that the above analysis is evidence 
for nature selecting the KMT solution, instead of the one 
presented here, but we wish to argue precisely the oppo- 
site. In the non-dynamical framework, there is nothing 
to suggest that the KMT solution is more valid or less 
valid than the new solution. The fact that two different 
observables are obtained is due to the immense freedom 
present in the choice of CS scalar, which points at an in- 
completeness of the non-dynamical framework. This is to 
be contrasted with the dynamical framework, where the 
coupling field is determined by its evolution equation and 
there is only freedom in specifying its initial condition. 



APPENDIX B: ARBITRARILY FAST-ROTATING 
BH SOLUTIONS 

In this appendix we study the possibility of extend- 
ing the analysis presented in this paper to a background 
that is not necessarily slowly-rotating. We shall see, how- 
ever, that dropping the slow-rotation approximation is 
sufficient to prevent us from being able to find analytic 
solutions. 

Consider then the Kerr metric in Boyer-Lindquist co- 
ordinates: 

,2 / 2Mr\ , 2 4Marsin 2 <9 ,, E l2 
ds 2 = - I 1 — ) dt 2 - dtdcf) + —dr 2 



+ Ed(9 2 +sin 2 



E 

2Mra 2 sin 2 9 



(Bl) 



where E = 



2Mr + a 2 . The 



perturbation shall be parameterized in the same way as 
before, leading to the following full metric: 



ds 2 



1 



2Mr\ , , , 

— J (l + h)dt 2 



AM a r sin 



(1 + uj)dtd(j> 



+ — (1 + to) dr 2 + E (1 + k) d9 2 
o 2Mra 2 sin 2 9 



(B2) 



All arbitrary metric perturbations, h, u>, to, k and p, are 
assumed to be at least of 0(C)- 

We can now study the evolution equation for the CS 
coupling field in the small-coupling approximation. In 
this limit, the evolution equation becomes 

Atf rr + 2 (r - M) i? >r + cot 9 t? >e + -d,ee 



24a 



a M 2 ) (rcosfl) 



(r 2 — 3a 2 cos 2 



E 5 



(B3) 



The general solution to any such inhomogeneous differ- 
ential equation will be given by the linear combination of 
a homogeneous and a particular solution. The homoge- 
neous solution is 



$H = #H(r)#H (0), 



(B4) 



where 



EiF 



E 2 F 



a a 
2' 2 

I I 
2' 2 



,ot,x 



,0,X 



{r-M + 1 )- &/2 
(r-M + 7)-^ 2 , 



0(0) = E 3 P_ &/2 (cos 9) +E 4 Q_ &/2 (cos 9), (B5) 
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where again P(-) are Legendre polynomials of the first 
kind, Q(-) are Legendre polynomials of the second kind, 
F(-) are generalized hypergeometric functions, Ei are 
constants of integration and the coefficients a and (3 are 
given in Eq. ([33]) . We recognize this solution as a gener- 
alization of Eq. (|3^|) to arbitrarily rotating BHs, where 

7 = y/M 2 - a 2 , 



Due to the strong similarities between this solution and 
the slow-rotation one, we directly conclude that the con- 
stants of integration must be chosen such that — 



const., by requiring asymptotic flatness, reality and a 
well-defined stress-energy tensor. 

The problem then reduces to finding the particular so- 
lution to the evolution equation for arbitrarily fast BH 
rotation. Unfortunately, such an equation does not ap- 
pear to be amenable to symbolic manipulation. Without 
such an analytic solution, we cannot proceed with the 
boot-strapping algorithm described in the main paper. 
However, the differential operator associated with the 
Ricci and C-tensors for the metric in Eq. (|B1[) is truly 
formidable. Thus, the probability of finding a solution 
for the metric perturbation, even given an analytic form 
for d, is small. 
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